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We explore the consequences of assuming that the neutrino mass matrix is a linear combination 
of the matrices of a three dimensional representation of the group S3 and that it has one zero 
mass eigenvalue. When implemented, these two assumptions allow us to express the transformation 
matrix relating the mass eigenstates to the flavor eigenstates in terms of a single parameter which 
we fit to the available data. 

PACS numbers: 14.60. Pq 



1. INTRODUCTION 

Following the discovery of neutrino oscillations, there has been considerable progress in determining values 
for the neutrino mass differences mf — m| and for the mixing angles relating the mass eigenstates to the flavor 
eigenstates. The most recent fits suggest that one of the mixing angles is approximately zero and another has 
a value that implies a mass eigenstate that is nearly an equal mixture of i'^ and i^r- If these conclusions were 
exact, then they could be accommodated by postulating a neutrino mass matrix having a symmetry based on 
a three dimensional representation of the permutation group S3. This connection has been extensively studied 
in the papers listed in Ref. [l|. The approach taken here is to retain a remnant of the S3 symmetry, assume 
that one neutrino mass is zero, and see what this implies about the final form of the neutrino mass matrix and 
transformation between mass and flavor eigenstates. 

For Majorana neutrinos the most general form of the mass matrix is 

A Bi B2 

AU ^ \ Bi Ci D \ (1) 

B2 D C2 

Experiment seems to show approximate fJ. — t symmetry in the sense that one mass eigenstate has an almost 
equal probability of being or Vr- To realize this with Mi, requires Bi « B2 and Ci « C2. As mentioned 
above, exact — t symmetry can be nicely modelled using a 3-dimensional representation of the finite group 
S3. However, suppose fJ- ^ t symmetry is not exact but we assume can still be expressed by the matrices of 
53. This ansatz, together with the assumption that one of the neutrino mass eigenvalues is zero, as required, for 
example, by the minimal seesaw model, allows us to derive two relations among the mixing angles and to predict 
all of the mixing angles in terms of one parameter. 

In the next section we review the conditions imposed by 6*3 on the elements of Al^,. In Sec. 3 we discuss the 
effect of these conditions on the minimal seesaw model. Following that, in Sec. 4, we find the eigenvalues and 
eigenstates when one mass eigenvalue is zero. Then, in Sec. 5, we write A, . . . , _D of M^, in terms of the mixing 
angles in the usual way and use the conditions derived from ^3 and from having one eigenvalue zero to find 
relations among the mixing angles. We are able to express these angles in terms of one parameter. In the last 
section we summarize our conditions on the mixing angles and compare our predictions with experiment. Finally, 
in an Appendix we discuss why it is possible to study neutrino mixing separately from the charged lepton sector. 
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2. CONDITIONS ON THE MASS MATRIX FROM S3 



The three dimensional representation of 6*3 is weU known. Nevertheless, for clarity, we will repeat it here. 
Each line of the following gives the elements that belong to a particular class 

/I 0\ 

D{e) =010 (2) 
\0 ij 

/00l\ /oio\ 
D(a) = 10 0, D{b) =001 (3) 






D{c) = 1, D{d) =010, D{f) =100 (4) 



This is a reducible representation and the sum of the elements in each class commutes with every element of the 
group. If we define 

Di = Die) , D2 = D{a) + D{b) , D3 = D{c) + Did) + D{f) (5) 
then the most general mass matrix, invariant under S'3, is 



M = aDi+l3D2+jD3=\BAB\ (6) 




Clearly this is not general enough but we can get a matrix which still respects ^ — t symmetry by breaking S'3 
with £)(c), 

M = aDi+pD2+jD{c) ^ \ f3 a P + j ] . (7) 

V /3 /3 + 7 « / 

This has the additional condition A + B — C + D necessary for tri-bimaximal mixing. 
But suppose n — T symmetry is not exact so we break S'3 in all possible ways 

/ a + -f P + e l3 + 5\ 

M = aDi+ liD2+-iD{c)+5D{d) + eD{f) = /3 + e a + 5 /3 + 7 (8) 

+ /3 + 7 a + ( j 

where we don't include D{a) or D{h) because we would have to add them to get a symmetric matrix and their 
sum is I?2, and we omit I?3 because it just adds the same amount to each matrix element. Thus we get something 
of the form of Eq. ([T]) but the important thing is that however we break the /i — r symmetry there are still two 
relations among the elements, 

2A + B1+B2 = C1+C2 + 2D (9) 
= C2-C1. (10) 

These remnants of the S3 symmetry are what we will use to restrict the parameters of the minimal seesaw model 
and to restrict the texture of the neutrino mass matrix. We will refer to them as the S3 conditions. 



3. MINIMAL SEESAW MODEL 

In the minimal seesaw model (as reviewed, for example, in Ref. 0]) the neutrino mass matrix is written as 

= mDM^^m'j^ (11) 
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where tou is a 3 x 2 matrix 



mo = I &i ci 

h C2 



and Mfj is a 2 x 2 matrix, 



Mr 



M22 M23 

M23 Af33 



(12) 



(13) 



The rank of M^, is two and therefore one of the eigenvalues of ([TT|) must be zero. 
Evaluating ([TT|) for the parameters in ([T|) gives 

^ a2M22 - 2aia2M23 + af M33 



2(^1+52) 
D 



i(Ci-C^2) 



a2(ci + C2)A'f22 - [a2('>i + ^2) + ai(ci + C2)]M23 + ai(6i + 62)^^33 

2P 

(cf + cl)Af22 ^ 2(6iCi + b2C2)M23 + (fci + bl)M33 

2V 

C1C2M22 - (&2CI + &lC2)M23 + &l62A^33 



V 

a2(ci - C2)Af22 - [02(^1 - ^2) + ai(ci - C2)]M23 + ai(6i - 62)^ 



33 



2P 

(cf - cl)Af22 - 2(6ici - b2C2)M23 + {bl - b|)M33 
21? 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



where V — M22M33 — A'f|3 is the determinant of ([T^ . Now using the 5*3 conditions © and (flUl) we get, after a 
lot of simplification. 







(2a2 - ci - C2) 

- (2ai-6i-52) 
= (C1-C2) 

- (&1-62) 



(02 + ci + C2)—; (ai +61+62) 



(02 + Ci + C2 



2? 

^M2a 
' V 



(ai + 61 + 62) 



V 
V 



M22 , , N Af23 

[0.2 + ci + C2)—; (ai + 61 + 62) 



(02 + ci + C2)^ (ai + 61 + 02)^- 



(20) 



(21) 



Since we assume that 61 7^ 62, c\ ^ C2 (to avoid Bi = B2, Ci ~ C2), Eqs. (|20p and ((2T|) could be solved by 
requiring 

ai + 61+62 = 0, (22) 

a2 + ci+C2 ^ 0. (23) 

But, when we put these relations back into (fH)) - (|19p . we get additional unwanted constraints 

B1+B2 = -A, (24) 
C1+C2 + 2D = A. (25) 

Thus the solution of ([20l) and ([2T|) must involve conditions on the parameters of A/fl as well as those oi mi). 
Another way to understand the S3 conditions and the restrictions ((22)) . ((23)) is to consider a Z2 symmetry 



ii 



(26) 
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where 



2 k 2 2 

= I 2fc fc2 _2 I . (27) 



2k -2 fc 
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Eq. (|26|) gives two conditions 



B1+B2 k 

(28) 



Ci + C2 + 2D - 2yl 

Bi - B2 1 



(29) 



which are the S-i conditions if fc = — 1 . 



In the fohowing sections we turn to finding restrictions on the mixing angles. In those sections the only use 
we will make of the minimal seesaw model is as motivation for setting one mass eigenvalue equal to zero. 



4. EIGENVALUES AND EIGENSTATES 

For a zero mass eigenvalue we have, 

B, Ci D /? = A /? = 0. (30) 
B2 D C2 ) \l ) \l I 

If we assume a ^ (we will check this below) then we get three equations 

A = -PB1-GB2 (31) 
Bi = -pC'i-aD (32) 
B2 = -pD-aC2 (33) 

where p = P/a, a = ^/a. So -B2 and A are given by ([32]) . ([33|). and 

A = p^Ci + a^C2 + 2pa D (34) 

Now let's use these in the S3 relations. Eq. ([TU]) and Eq. © give 

{a-l)C2-{p-l)Ci + ip-a)D = (35) 
(2p2 - p - l)Ci + (2(7^ -a- l)C2 + (4pcr - p - a -2)D = Q (36) 



Eqs. (|55|) and can be reduced to 



If we choose the solutions 



{p + a + - cj)C2 + [l - p)D] = 0, (37) 
{p + cT + l)[{l- p)Ci + {l-n)D] = 0, (38) 



Ci = -^D, (39) 
1 - p 

C2 = -^D, (40) 



then — C1C2. The mass eigenvalues that we expect to be nonzero are given by 

1 



m± 



A + Ci+C2± V{A + Ci+ C2Y + 4(p2 + ^2 + 1)(Z?2 _ C1C2) (41) 
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where we have used ([5^ and Thus this solution makes a second mass eigenvalue zero. We need nonzero 

two masses in order to have two oscillation lengths. We might tolerate two zero masses in the case of normal 
hierarchy, ^ rn2 ~ rni. We will ignore that special case except for a brief comment at the end of Sec. 5. 
Thus the only way to avoid two zero masses is to require p + a + 1 = 0. 

If we take Ci and C2 as the independent variables the nonzero eigenvalues are, from (|4T|) using (|34p . 



^± = 2(2^) ^2 



(2 + 2(7 + (t2)Ci - (l + a2)C2 



2o-+ 1 



The total set of eigenvalues and eigenfunctions can be reduced to 

mo = 0, 



Wo > 



1 



2^1 + Re{a) + \a\^ 
1 



\Ve > -(1 +Cr)|i^^ > +<T\Vr >] 



nia 



2(7+1 



[{a + 2Y Ci~{cr- 1)' C2 



1 



Wb > 



2(^2 + 1) 



2ct + 1 



[C2 - Ci] 



1 



V6^/l + Re{a) + \a\ 



-(1 + 2a)|i/e > -(1 - > +(2 + a)Wr >] 



where a, b are +, — or — , + depending on the signs of the factors in the absolute value of (|42p . 
The only remaining case is to go back to (pO| and set a — 0. Eqs. (|31l) . (|32t . and (|33p are then 



where A = /3/7. The ^3 conditions are 



B2 
D 
C2 



Bi 
A 



-XBi 
-XCi 

-XD = X^Ci 



(A - l)Ci 



Since one eigenvalue is zero the remaining eigenvalues are given by 



Ad 



A + Ci+C2±JiA + Ci+ C2f + 4(i?2 + Bl + Bl - C1C2 - ACi ~ AC2) 



(42) 

(43) 
(44) 

(45) 
(46) 

(47) 
(48) 



(49) 
(50) 
(51) 



(52) 
(53) 



(54) 



and dm) - ([53]) give D'^ + BI+B2-C1C2- ACi - AC2 = . Thus a = would require a second mass eigenvalue 
to be zero. 



So the only solution with broken p — t symmetry and two nonzero masses is given by (|43p - ()48|) . Since we 
know fJ- — T symmetry is approximately true the parameter a will need to be large for inverted hierarchy or 
approximately — i for normal hierarchy. 



5. RESTRICTIONS ON THE MIXING ANGLES 



The conditions on the elements of the mass matrix A, . . . ,D will allow us to put conditions on the mixing 
angles {9s, da, 6x) = {O12, 823, dis)- The neutrino mixing matrix Q which diagonalizes M^, via V^Mi,V = M^'^s 
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can be decomposed as V = U"UU' Q where C/ is a CKM type matrix 



U — I SsC-a CsSaSxP CsCa ~\~ SsSaSxP 




(55) 



with (scjCq) = (sin^Q, cos^) for a = s,a,x and p = e"^^" where Sjj is the Dirac phase, U" is the rephasing 
matrix, diag(e'"i , e*"^ e*"3)g, and U' = dia.g{e~"'''/^,e-''f'^/'^ ,e'"^^^^) where 0i, (j}2, and 03 are Majorana 
phases. The neutrino mass matrix is then 



with elements given by 



A = 


[ 2 2 1 , 2 2 1 , *2 2 M -2iai 




Bi = 


\px [-^fiCgCa pSxSaC^s\^^l \^s^s^a ~l" P^x^a^ s 


m'2 +p*CxSxSam'j] e 


B2 = 


\px [-^fiCgSa -j- pSxCaCg^ni^ Cx [SsCsSa P^x^a^ s 


m'2 ~ p* SxCxCam'.^] e 


Ci = 


[{SsCa - PSxCsSafm^ + {CgCa + pS^SsSafm^ 4- 






[{SgSa +pSxCsCafm[ + [CgSa - pS^S sCa)"^ m'2 -+ 




D = 


[{SsSa + pSxCsCa){SsCa - pSxCsSa)m[ 





-i(ai+Q2) 



(56) 

(57) 
(58) 
(59) 
(60) 
(61) 



CsSa - pSxSsCa){CsCa + pSxSsSa)m'2 - dsaCam'^e 



-1(02+03) 



(62) 



where 
moVT" 



moVT+Te^'^i, m'2 = moe*' 



for inverted hierarchy or m'l = 0, m'2 = mo-s/re"^^, 77x3 



'''^ for normal hierarchy. This mo is a universal mass, not the same as the eigenvalue of the previous 
section, and r is the ratio of the mass splittings, r = A^/Aa. It is easy to see that n — t symmetry requires 
Sa = Ca and Sx = 0. Tri-bimaximal symmetry requires, in addition, tan^s = \/2 or — l/\/2- 

We have four relations for a from pi]) . (15^ . ([55]) . and (|55|) . all with p replaced by —a — 1. The first three are 



(63) 
(64) 
(65) 



1 


Bi 


-B2 


a 


A 


-Bi 


1 


Ci 


- D 


a 


Bi 


-Ci 


1 


D- 


~C2 


a 


B2 


- D ■ 



Now we assume inverted hierarchy and substitute (|57p - (|62p on the RHSs to find relations on the mixing angles. 
Only two of these are independent because we have set rrig equal to zero so we get just one relation among the 
angles. 



S T" C — Cx 



g gi(a3-ai) _ g gi(Q2-Ql) 



(66) 



and the solution for cr, 



1 

(7 



„i(Q2 — 03) 



(67) 



One of the mass eigenvalues m^, mf,, given by (145^ and (j47p . must equal m'^ and the other mj. If we evaluate 
mo and mb using (|67| and Ci, C2 given by (|60| . (f6T|) we find this requires a2 = as. Thus cr is real and (|66| 
requires 



(5^1 + — Q!2 = or TT 



(68) 
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It is immaterial which sign we take from (|68p: in what follows we use 



Sx = Cx{Ca-Sa) (69) 
1 ^ ^^Q^ 



We still have the condition from 5*3, Eq. which now depends only on Sd, 

1 _ 2L> - Ci - C2 
a ^ 2C1-C2-D' 

Using ([70)1 for the LHS and ([66| for e'''° , this gives a quadratic equation for tan 9s 

mie — 7712 e 



(71) 



tan0. 



(sa + Ca)(rn1 + — 2miTO2 cos 2a) 



Ca;(l - ACaSa){mi - 7712) ± \J cl{l ~ ACaSaY{mi - 1712^ +4(1 + 2 Cq Sa)("T-i + m| - 2miTO2 COS 2a) (72) 

where a = ai — Q!2- This gives real values only if 5]j is zero or tt. One solution of the quadratic equation is 

tan 6*8 = — ^[I-SqCo] (73) 

Ca + Sa 



or, when we use 



The other solution is 



tan^, = ^ (74) 



tan 6*5 = [l+2CaSQ] = C^(Ca + Sa). (75) 

Since tan(-| — Og) — l/tan^s, and since the oscillation experiments measure sin^ 29, these two solutions are 
effectively equivalent. 

Another way of expressing the results is to write all of the mixing angles in terms of the one parameter a. 
From ^ and then ^ we find 

tan 0a = (76) 
a 

tan 6*^ = = , (77) 

Vl + 2(7 + 2ct2 ^ ^ 



and, from ([73]) . 



or, from (j75p . 



^/2V1 + g + a2 

tan 61, === — — (78) 

i + Z(T 



tan0, = ^^l±i^=. (79) 



These two solutions have — m[ — movT+T, ma — m'2 — mo for (|78p or = = mo, ma = m[ = 
moy/1 + r for ([75)1 . (Here we are only concerned with the magnitude of the masses so we have neglected 
Majorana phases.) Since \a\ is large, Eq. (|48l) shows that \i^b > has a larger fraction of l^e > than does > so 
mb should be smaller than ma- If r > then Eq. (|75p gives ma ~ moVT+r > mt = mo but Eq. (|74p would 
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require r < 0. However, as mentioned above, these two solutions for tan 0s are effectively indistinguishable so 
we won't worry about this point further. 

So far we have considered only inverted hierarchy but Eq. (j44p supports normal hierarchy with cr « — If we 
set mi = and compare the right hand sides of (1551) . ([M)) . and ([55]) we get the condition 

CsiCx - Sx (Sa - Ca)) + Ss(Sa + C^) = , (80) 

where to make the algebra simplier we immediately neglect all the phases. This expression for tan^g, 

tang, ^ ^ C.-gx(5a-Ca) ^ ^g^^ 
~t~ Ca 

superfically looks different than those of the inverted case. However, if we now set d7T]) equal to or (IMl) or 
(l65t we find Eq. (l69t which, when combined with (IHlT) , reproduces our first solution for tan 0s in the inverted 
hierarchy case, Eq. (f74| . 

Using these conditions the solution for a is now 



2s„ 



or 



(82) 



tanOa - J-^ (83) 



When we solve for nia, mb of Eqs. P5)) . (|T7)) we get 



TOft = TO3 = moVT+r (84) 

JTla = "^2 = TO0\/?' (85) 



If we set TO2 = rather than mi, we get the second expression for tandg of inverted hierarchy, (|75p. and ma = mi. 
For this normal hierarchy case a is approximately — ^ and the eigenf unction with zero mass, Eq. (|44|) . has the 
largest fraction oi \ ve > as it should. 

We have insisted that only one mass eigenvalue be zero. If two masses are zero then, from ([57l) - (|62|) with 
only mg nonzero, the 6*3 conditions are both satisfied by (j66|) above. 



6. SUMMARY 



We assume that the neutrino mass matrix is given by the three dimensional representation of the group S3 and 
that one (but only one !) of the mass eigenvalues is zero, as required, for example, by the minimal seesaw model. 
The experimentally testable results for either inverted or normal hierarchy a real Dirac phase factor (there could 
be nonzero Majorana phases but they don't affect our results), 

6d = ,7r (86) 

and two conditions on the mixing angles, 

tan 9x = Ca - Sa , (87) 

and 

tanOs ^ (88) 

Cx [Sa I Ca ) 

or 

tan0s = Cx{Sa + Ca) . (89) 
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For inverted hierarchy we can convert (175)) - ([75]) to 



and 



or 



sin^g.^ f + " + (92) 

3(1 + 2cr + 2cr2) ^ 

sin^e, = _ii±^^ll_. (93) 
3(1 + 2c7 + 2ct2) ^ ' 

to more easily compare with the data. The result of fitting these expressions to the experimental values 0, [1] 
is shown in the following table. As mentioned above the oscillation expressions depend on sin^ 29 and thus the 
experiments can't distinguish between 9s greater or less than 7r/4. So the fit using (j92p assumes the experimental 
value is less than 7r/4, while that using (|93l) assumes it is greater than 7r/4, since these are the values indicated 
by the formula. (Approximate /i — t symmetry gives Sa ^ Ca ^ "71 (|88p gives tan 9s ^ while (I89|) gives 

tan^s \/2; thus a fit of (|93|) with 9s < 7r/4 is untenable.) The two sets of expressions give equivalent fits to 
the data as they must. 



Angles 


Best Fit 


Exp Range 


Fit 1 


Fit 2 


sin^ 6a 


0.466 


0.408 - 0.539 


0.425 


0.425 


slv? Ox 


0.016 


0.006 - 0.026 


0.011 


0.011 


sin^ 6s 


0.312 


0.294 - 0.331 


0.337 




sin^ 6s 


0.688 


0.669 - 0.706 




0.663 



TABLE I: The second column gives the experimental best fit, the third column gives the 1 — a experimental range, the 
fourth column gives the central values using (|92ll . the fifth column gives the central values using (|93|l . The minimum 
is 2.15. The value of a which gives the minimum is —7.17 if we use (|70|) or —0.388 if we use H82[). Our fits give values for 
the angles of \ea\ = 40.7°, \6x\ = 6.02°, and |6I,| = 35.5° or 54.5°. 



Similarly for a normal hierarchy of masses we could find expressions for sin^ 9a, sin^ 9^, and sin^ 9s in terms 
of the a for normal hierarchy given by ([5^ . But if we call that ctat, and the a given by (|7D|) aj, then ([7^ and 
(IMl) give 

™ = 

and if we replaced ctat by c/ in normal hierarchy expressions for sin^ 9a, sin^ 9x, sin^ 9s we would reproduce (|90)) 
- (|93l) . The normal hierarchy expressions would be just a reparameterization of (|90|) - (|93|) above and thus give 
an identical fit. 

Recently MINOSS has presented a measurement of sm^{29a) sm^{29j:). We can easily fit this assuming the 
values above for sin 9a and sin^ 9s but replacing sin^ 9x by this combination. This is shown in Table II. The fact 
that the fit values of sin^ 9x are smaller than the value in Table I despite the experimental number being bigger 
is because of the larger error in that number. Still as increases by one from its minimum, sin^^; varies from 
0. to only 0.024, which implies that this model prefers small 9x. 
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Angles 


Best Fit 


Exp Range 


IH Fit 


NH Fit 


sm Ba 


0.466 


0.408 - 0.539 


0.449 


0.443 


sm2(26'a)sin2(26»,) 


0.18 


0.06 - 0.32 


0.021 




sm2(26»a)sm2(26»a;) 


0.11 


0.04 - 0.21 




0.025 


sin^ Qs 


0.312 


0.294 - 0.331 


0.335 


0.335 


Minimum 






2.90 


2.37 



TABLE II: The fourth and fifth columns give the results of fitting the MINOS values given in row two or row three. Again 
the experimental range is 1 — cr. The MINOS numbers depend on whether they assume inverted or normal hierarchy. 
The fitted numbers correspond to a sin^ of 0.0052 for IH or 0.0064 for NH. The values of the angles are therefore 
\ea\ = 42.1° or 41.7°, \e^ = 4.14° or 4.59°, and |6'J = 35.4°. 
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Appendix A: Charged Lepton and Neutrino Sectors 



In recent years it has become common to attempt a unified treatment of the charged lepton sector and the 
neutrino sector. In this paper, we discuss only the neutrino sector, as do many of our references. This appendix 
shows that the sectors can be discussed separately. 

The essential point of fiavor mixing is that the physical mixing matrix, being misaligned between two repre- 
sentations of the up-type and down-type fermions, is independent of formalism or representation. For the lepton 
sector, the physical mixing matrix is the so-called PMNS matrix (5|, 



PMNS 



(Al) 



where the two flavor mixing matrices are determined by 

UlM.MlU, = D,Dl , 



(A2) 



As in the body of the paper, we take the neutrinos to be Majorana particles. The diagonal mass matrices are 
denoted as Df, and D^, for charged leptons and neutrinos. Now, we can make an arbitrary rotation on all the 
lepton fields, including left-handed charged leptons and neutrinos as well as the right-handed charged leptons. 
Since the left-handed charged leptons and neutrinos reside in common SU(2)l doublets, they share a common 
rotation. 



ink 



(A3) 



Then the charged lepton and neutrino mass matrices become 



^ M, = TlM.Tl , 



(A4) 



and we denote the modified mixing matrices as Ue and Ui, respectively. Eq. (jA2l) becomes 

UImMu, = D,Dl , UjM.U, - , 



(A5) 
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where 

Ue = , = TlU, . (A6) 

The important point is that the physical mixing matrix is not affected, 

VPMNS = UlU^ = uItITlU, = UlUu = VPMNS ■ (A7) 

This is expected because, if it were not true, the physical mixing matrix would depend on the formalism or 
representation. This property of formalism/representation independence allows us to rotate the charged leptons 
to a mass diagonal basis, since what we want to discuss is just the physical mixing matrix. It doesn't matter in 
which basis the discussion is made. The question is how to realize this. 

We should also note that, after gauge symmetry breaking where the fermions acquire mass, the up- type and 
down-type fermions' mass matrices should be constrained by different representations of some symmetry if there 
is any. Otherwise, the two mass matrices would be constrained to be of the same form and this would lead to 
trivial physical mixing. The full group is at least a product. 

We can imagine that the 6*3 discussed in current work is kind of residual property of some symmetry. Before 
symmetry breaking, there would be a larger group governing both the charged lepton and neutrino sectors, 
especially the left-handed ones since they share a common left-handed doublet. But experimentally the symmetry 
is broken. There can be a residual symmetry for the left-handed charged leptons, for example -Ea, 



(1 

Z3 = {I,F,F'^} with F = 

V 



, ,2 



(A8) 



where w = e^"/^. If we use to denote the group elements of Z3, then the charged lepton's mass matrix has to 
satisfy 

F^M^MIF = DeDl . (A9) 

It can be verified that under this constraint MeMj has to be diagonal, 

M,Ml = D,Dl , (AlO) 

and Ue ~ I ■ In other words, the physical mixing comes solely from the neutrino sector. Actually, Z3 is a 
subgroup of 53 so we can apply ^3 in the charged lepton sector. But, as argued above, the representations of the 
charged lepton and neutrino sectors should be different. In other words, the residual Z3 of the charged lepton 
sector cannot be simply embodied in the residual 6*3 of the neutrino sector. The unified group should be at least 
a product group [l3| 

g = z^®s^. (All) 

If we want to apply 6*3 in the charged lepton sector too, this can be achieved by embedding in another 6*3 
whose representation is different from that of neutrino sector. Then the product group would be 3^(^)3^. 



[1] S. Pakvasa and H. Sugawara, Phys. Lett. B73, 61 (1978); 82, 105 (1979); E. Durman and H.S.Tsao, Phys. Rev. 
D20, 1207 (1979); Y. Yamanaka, H. Sugawara, and S. Pakvasa, Phys. Rev. D25, 1895 (1982); K. Kang, J. E. 
Kim, and P. Ko, Z. Phys. CT2, 671 (1996), hep-ph/9503346 ; K. Kang, S. K. Kang, J. E. Kim, and P. Ko, Phys. 
Lett. A12, 1175 (1996), hep-ph/9611396, M. Fukugita, M. Tanimoto, and T. Yanagida, Phys. Rev. D57, 44299 
(1998), hep-ph/9709388; H. Fritzsch and Z-z Xing, Phy s. Rev . D61, 073016 (2000), hep-ph/9909304 , E. Ma and G. 
Rajasekaran, Phys. Rev. D64, 113012 (2001), |hep-ph/010629l| P. F. Harrison and W. G. Scott, Phys. Lett. B557, 76 



12 



(2003) arXiv:hep-ph/0302025l ; S.-L. Chen, M. Frigerio, and E. Ma, Phys. Rev. D70, 073008 (2004) [Erratum-ibid. 
D70, 079905 (2004)] hep-ph/o404084; F. Caravaglios and S. Morisi, ^-ph/0503234; W. Grimus and L. Lavoura, 
JHEP 0508, 013 (2005), hep-ph/Q5Q4153: J. E. Kim and J. -C. Park, JHEP 0605, 017 (2006), hep-ph/0512130[ 
R. N. Mohapatra, S. Nasri, and H. B. Yu, Phys. Lett. B639, 318 (2006), hep-ph/0605020: R. Jora, S. Nasri, and 
J. Schechter, Int. J. Mod. Phys. A21, 5875 (2006), hep-ph/0605069 ; M. Picariello, Int. J. Mod. Phys. A23, 4435 
(2008), hep-ph/0611189; Y. Koide, Eur. Phys. J. C50, 809 (2007), hep-ph/0612058 , A. Mondragon, M. Mondragon, 
and E. Peinado, Phys. Rev. D76, 076003 (2007), arXiv:0706.0354 (hep-ph); A. Mondragon, M. Mondragon, and E. 
Peinado, AIP Conf. Proc. 1026: 164 (2008), , arXiv:0712.248 8 (hep-ph); C.-Y. Chen and L. Wolfenstein, Phys. Rev. 
D77, 093009 (2008), arXiv:0709. 3767l 

[2] W-1. Guo, Z-z. Xing, and S. Zhou, Int.J.Mod.Phys. E16, 1 (2007). 

[3] S.-F. Ge, H.-J. He, and F.-R. Yin, arXiv: 1001. 0940 

[4] S.-F. Ge, D.A. Dicus and H.-J. He, in preparation. 

[5] B.Pontecorvo, Sov.Phys.JETP 6, 429 (1958); Z.Maki, M.Nakagawa, S.Sakata, Prog. Theor.Phys.28, 870 (1962). 

[6] A. Barroso and J. Maalampi, Phys. Lett. B132, 355 (1983). 

[7] G.L.Fogh, E. Lisi, A. Marrone, A. Palazzo, A. M. Rotunno, Phys. Rev. Lett. 101, 141801 (2008) |arXiv:0806."2649] 

and arXiv:0809. 2936 hep-ph]. 
[8] G.L. Fogh et. al, arXiv:0805.25l7v3 [hep-ph] and Phys. Rev. D78, 033010 (2008). 

[9] P. Adamson et. al, Phys. Rev. Lett. 103, 261802-1 (2009). 

[10] C. S. Lam, "The Unique Horizontal Symmetry of Leptons", Phys. Rev D 78, 073015 (2008), ' arXiv:0809.T785] [hep- 
ph]]. 



